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1. Introduction 

Optical fiber grating sensors are becoming increasingly 
widespread in a wide range of applications^. Their char- 
acteristics include immunity to electromagnetic interfer- 
ence, high sensitivity, fast response, and keeping their 
calibration with time. Furthermore they possess a high 
tolerance to harsh environments such as sea water, con- 
crete, and extreme temperatures. The measured parame- 
ter, the resonant wavelength of the grating, is usually lin- 
ear in the magnitude of the measurand, e.g. temperature 
or strain, allowing for an accurate quantitative estimate 
of the measurand. 

Photonic crystal fibers (PCFs}^ are optical fibers, in 
which the cladding consists of a microstructured array 
of air holes running along the fiber axis. Here we con- 
sider index guiding PCFs, which have a solid core in the 
form of a missing air hole in the center. In index guid- 
ing PCFs light is confined to the core by modified total 
internal reflection, analogous to guiding in standard opti- 
cal fibers. The holes in the cladding are most commonly 
placed in a periodic triangular structure, which is char- 
acterized by a lattice constant (pitch) and the diameter 
of the holes. PCFs are typically made of silica^ or, more 
recently, polymer materials, where they are termed mi- 
crostructured polymer optical fibers (mPOFs)^. The hole 
structure in the cladding determines the optical proper- 
ties of the PCF, which allows for a large degree of freedom 
in tayloring key properties, such as dispersion and group 
velocity. 

In standard optical fibers the silica in the core contains 
dopants to increase the refractive index. The refractive 
index of the doped silica is sensitive to UV radiation and 
this may be used for inscription of gratings. This tech- 
nique can not be used for pure silica PCFs, since the core 
is not UV sensitive. For long-period gratings (LPGs) in 
PCFs other techniques exist, such as using an electric 
aroi, a CO2 laser—, or mechanical microdeformation^. 
Otherwise the PCF can be specially made with a UV 
sensitive core for fiber gratings allowing for the inscrip- 
tion of both LPGs and Bragg gratings (BGs)i. 



Polymer has a smaller Young's modulus and greater 
elastic limit than silica, which makes it a good material 
for measuring strain and bending. In addition, polymer 
is biocompatible and can bind biomolecules directly to its 
surface, making it ideal for biosensing. Standard poly- 
mer optical fibers are made of poly-methyl methacrylate 
(PMMA), which is inherently photosensitive, allowing for 
both bos'* and LPGs^ to be inscribed using a UV lamp. 
The inherent UV sensitivity was also recently used to 
write BGs in mPOF^*', whereas heat treatment under 
mechanical stress was used to write the first permanent 
LPG in an mPOF^i. 

An important property of grating sensors in silica 
PCFs are that they are largely temperature insensitive, 
showing - 1 pm/K for both BGsM and LPGsH. The 
strain sensitivity is of the same order for LPGs and BGs 
in PCFs, and thus strain sensors may be realized with 
reduced cross sensitivity to temperature sensing. In ad- 
dition, the holes of a PCF can be infiltrated with a sub- 
stance having a high temperature coefficient to reduce or 
enhance the temperature sensitivityi^. 

Whereas fiber-optic strain and temperature sensors use 
the material properties of the fiber, fiber-optic biosensors 
and refractive index sensors use the principle of evanes- 
cent wave sensing. In particular, there is a growing in- 
terest in optical biochemical sensorsi^. Here PCFs have 
the advantage that the analyte may be infiltrated into the 
air holes and thereby create a strong interaction with the 
probing electromagnetic fieldi^iii. When using mPOFs 
one can take advantage of different polymers with special 
biosensing properties, such as TOPAS^^. 

Liquid infiltrated PCF-LPCsi^'^o and PCF-BGsiii^ 
have been demonstrated, which is an important step to- 
wards the ultimate goal of developing a so-called label-free 
biosensor in a PCF. Fiber gratings potentially allow for 
label-free biosensing with PCFs^S, by tracking how the 
grating resonance changes during operation. This may 
open for new possibilities in biomedical applications^. 

The dispersion plays a crucial role for all types of grat- 
ing sensors in standard optical fibers^^i^l. The dispersion 
of PCFs depends strongly on the hole structure in the 
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cladding and thus the performance of the PCF grating 
sensor can be designed by an appropriate choice of pitch 
and hole size. 

In this paper we explore the possibility of optimizing 
the hole structm'c of PCFs for sensing purposes. We 
do this by analyzing the sensitivity of triangular PCF 
gratings for different pitch and hole diameters. The ap- 
plications under study are biosensing, refractive index 
sensing, strain, and temperature sensing. We show that 
the sensitivity is strongly influenced by the choice of the 
structural parameters of the PCF. We employ realistic 
values for the fibers parameters and give realistic values 
of pitch and hole diameter that minimize or maximize 
the sensitivity to different measurands. 

We find that BGs are best characterized in the usual 
way by their sensitivity, 7, defined as the resonant wave- 
length shift divided by the resonant wavelength. In con- 
trast, we find that the LPGs are best characterized by 
their quality factor, Q, defined as the resonant wave- 
length shift divided by the full-width at half maximum 
(FWHM) of the resonance dip. 

2. Photonic crystal fibers 

The eigenmodes of the PCFs for a temporally harmonic 
electrical field, E(r,t) = Ei^(r)e~*"*, may be found from 
Helmholtz eigenvalue equation, 

^2 

V X V X E^(ri, z) = ^£(ri)E^(ri, z), (1) 

where uj is the angular frequency, c is the speed of light 
in vacuum, e is the dielectric function, and r^ — {x,y) 
is the transverse coordinate vector. The z-axis is chosen 
to coincide with the fiber axis. We look for plane-wave 
solutions of the form E^{r) = E^(rj_)e''^^ with /3 being 
the propagation constant. The effective index is defined 
as nefi(fc) = f3{k)/k, where k is the free-space wavevector 
and is related to the frequency and wavelength by fc = 
(jj/c = 2tt/X. The effective indices of the core mode and a 
resonant cladding mode are denoted by rico and rtci- The 
orthonormalization for Eq. ([T]) is 

/ dr^Ej„(r^)£(r^)Et (r^) = Smn- (2) 
Jn 

The PCF is typically made of two materials with a high 
refractive index contrast. A general dielectric function, 
e(r), for a two material triangular PCF structure is de- 
scribed by high and low index regions as seen in Fig. [TJ 
The periodic triangular structure is spanned by the lat- 
tice vectors i?i = (%/3, l)/2 and R2 = (%/3, -l)/2. The 
general dielectric function equals for the base region 
and for the hole regions, with nf, and being the 
refractive indices of the base and hole material, respec- 
tively. 

The effective index of the core depends strongly on the 
wavelength, the hole diameter, and the pitch. The hole 
diameter, d, and pitch, A, are indicated in Fig. [TJ Large 



hole diameters and small pitch also give a small core and 
therefore a strong effective nonlinearity. 

Triangular silica PCFs are single mode for d/A< 0.42 
for all wavelengths (25ll26j . For small hole sizes {d/A < 
0.3) the PCF is only weakly guiding, and bend losses 
increase drastically. Bend losses also become pronounced 
in the short wavelength limit, A ^ A, regardless of hole 
diameter. PCFs can be designed to have a large mode 
area. 

The transmission spectrum of an LPG in a PCF only 
shows a few cladding mode resonances, whereas several 
resonances are seen in a standard optical fiber LPG spec- 
trum. The cladding modes of a PCF are different from 
the cladding modes of a standard optical fiber, since the 
fiber cross sections obey very different symmetries. In 
particular a PCF obeys a six-fold rotational symmetry 
(C|) while the standard optical fiber has circular sym- 
metry. The effective indices of the cladding modes of a 
PCF are closely spaced, and this spacing is much smaller 
than the difference between the core effective index and 
the effective index of any cladding mode. In other words, 
the beat length between any two cladding modes is much 
longer than the beat length between a cladding mode and 
the core mode. In practice we can assume rid — '^fsm for 
all cladding modes, where the fundamental space filling 
mode, npsM, is the effective index of an infinite cladding 
structure spanned by the lattice vectors in Fig. [T] 

We have calculated the effective index for the core 
mode, Tico, and the fundamental space filling mode, 
?^FSM, for different hole to pitch and wavelength to pitch 
ratios, with the base material having a refractive index 
1.45, which is the refractive index of silica at the wave- 
length A ~ 1/im. The effective indices are shown in Fig [2] 
The effective indices decrease with increasing hole size 
(d/A) or increasing normalized wavelength (A/A). The 
effective index of the core mode is larger than that of the 
FSM, regardless of the fiber structure parameters (d,A). 




Fig. 1. A PCF with triangular cladding structure. The 
structure is characterized by the lattice vectors, Ri and 
i?2, the hole size, d, and the pitch, A. The refractive in- 
dices n/i and rif, are the hole and base material refractive 
indices 
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3. Coupled mode theory 



be expressed in a compact matrix form as 



A . Introduction 

Coupled mode theory^l is a simple and compact theory 
that has been successful in accounting for the character- 
istics of fiber gratings. It allows us to derive the quali- 
tative behavior of the gratings with some general, sim- 
plifying assumptions on the physical mechanisms. The 
theory assumes that near resonance the electromagnetic 
field propagating through the fiber may be expressed as 
a superposition of eigenmodes of the unperturbed fiber 
aco(z)Eco(r_L)e*'''^°^ -|-ai(0)Ei(r_L)e''''^, where 'z' denotes 
a resonant mode, and the a's are amplitude coefficients 
that depend on z due to the perturbation. For Bragg 
gratings the resonant mode is the counter propagating 
core mode, and for long-period gratings the resonant 
mode is a copropagating higher order mode, usually a 
cladding mode. The coupled-mode equations constitute 
a set of linear coupled differential equations. They may 
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Fig. 2. Effective indices for the core (top) and funda- 
mental space filling mode (FSM) modes for triangular 
photonic crystal fibers. The lines indicate different val- 
ues of the hole diameter relative to the pitch. The lines 
correspond to different values of d/A, increasing by 0.1 
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where the detuning in the resonance between the core 
mode and the resonant mode, j, is 



/3co - A - ^ 



27r 



and the coupling constant is 



= / dr^Etjrx)A£(r^)E,(r^), 
1 Jn 



(4) 



(5) 



where the integration is over the fiber cross sec- 
tion, 17. A£(r^) is the first phase matched 
Fourier component (±1) in the expansion Ae(rj^,z) = 
Em=-oo A£„(ri) exp(-im||z). 

In Eq. ^ we have neglected the diagonal terms, the so 
called self-coupling terms. The self-coupling arises from 
the perturbation of the individual modes due to the in- 
dex change of the grating. The self-coupling shifts the 
resonant wavelength some nm's, while otherwise leaving 
the spectrum of the grating unaltered. Since this effect 
is typically small we neglect it for our purposes. 

Using Eq. ^ we obtain a resonance condition in terms 
of effective indices and the wavelength 



A. 



(nco(Ar) - 

nf(Ai.)AG, 



.(Ar))Ac 



(6) 



where Ui is the effective index of the resonant mode. The 
equation is valid both for LPGs, fif = nco^nd, and BGs, 
Uf = 2nco. 

In the following we use the bar for indicating a differ- 
ence between the core and the resonant modes as done in 
Eq. ([6]). Thus ^ = f3co — Pi for the propagation constant, 
and for the group index we have ^ 



-\^dx{l/LB), where Lb 
, . = Ml 

'g.* — dk ■ 



27r//3 is the beat length, with 
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Fig. 3. Definition of the full- width half maximum 
(FWHM) in detuning, (5fwhm ^ 5^ — 5^ , and in wave- 
length, Afwhm = AA+ — AA^ . 
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The resonance conditions must be solved considering 
the dispersion of the effective indices, since these are 
strongly dependent on the pitch and hole size of the PCF 
structure as seen in Fig. [21 

B. Dispersion effects on the wavelength FWHM 

Couple-mode theory in the form of Eq. ([3]) will give an 
expression for the transmission in terms of the detun- 
ing and coupling constant. However, experimentally the 
transmission is measured as function of wavelength. 

In Fig.Othe FWHM for the detuning, Jfwhm = |<5+ - 
5^\, and the wavelength, Afwhm = [AA^ — AA~|, are 




Fig. 4. Beat length as function of the wavelength over 
pitch, A/A, for an LPG. The gray area indicates negative 
group index mismatch ng < 
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Fig. 5. FWHM as function of the wavelength over pitch, 
A/A for an LPG 



defined. By definition we have = T{5 ) = 0.5, and 

equivalcntly for the wavelength we have T{\ + AA+) = 
T(A,. + AA-) = 0.5. 

To find a relation between (5fwhm and Afwhm that 
incorporates the dispersion of the effective indices we ex- 
pand the detuning around the resonance wavelength to 
first order in AA 



5{K + AA) 



27r 

/3eo(A, + AA)-/3,(A, + AA)- — 



dkdX 



AA = -,lg(A,) JaA, (7) 



With this equation it is possible to obtain two equations 
relating and i5~ to AA"*" and AA~. Subtracting the 
two equations and taking the absolute value of both sides, 
we obtain an expression for (5fwhm and Afwhm 
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FWHM 



Ar 



27r |rig(Ar)| ' 



(8) 



which incorporates the dispersion of the effective indices 
in a consistent manner. 

We note that the expression is not valid at group index 
matching between the two resonant modes, fig = 0. In 
this limit, however, the FWHM goes to infinity, and any 
wavelength shifts of the sensor will become difficult to 
detect. 

C. Long-period gratings 

Long-period gratings couple the core mode to a co- 
propagating high-order mode. The difference in effec- 
tive index of two such modes is typically small, hence by 
Eq. ^ the grating period of such a grating must be long, 
Ag ~ 100 - 1000 X A. For LPGs we put the subscript i 
= FSM, and at = ofsm, <5 = Si = SpsM, k ^ Hi = kfsm 

Typically the FWHM of an LPG is of the order of 
nanometers. The minimal detectable shift in resonance 
condition is thus set by the FWHM and not the spectral 
resolution of the spectrum analyzer. 

The coupled-mode equations for the long-period grat- 
ing may be solved analytically from Eq. ([3]) yielding the 
transmission coefficient 



m = 1 - 



1 



1 



4|k|- 



■ sm 



(l^^lL^l 



4kF 



(9) 



At resonance (S = 0) the transmission is sinusoidal, 
T{6 = 0) = cos^(|k|L), where L is the length of the 
grating. The first zero for T is at \k\L — 7r/2. Ignor- 
ing the wavelength dependence of the coupling constant, 
K — K(Ar) = 7r/ (2i), we substitute ^ — S/{2\k\) and solve 
r(^) = 0.5 numerically This yields ^+ = ~ 0.80. 
We can find the FWHM of the detuning as (5fwhm/|k| = 
2(C+ - O ^ 3.2. The wavelength FWHM of the trans- 
mission can then be found using Eq. ([5]) 
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ng{X,)\ L 



(10) 
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The wavelength FWHM obviously is large when the 
group index matching is close to zero, «g(Ai-) = 0. 

In conclusion we have found an expression, characteris- 
tic for LPGs, that relates Afwhm directly with the length 
of the grating, the group index and the resonant wave- 
length. The expression is derived under the assumption, 
that \k{X)L = n/2 

D. Bragg gratings 

Bragg gratings couple the incident core mode to a 
counter-propagating core mode, hence (3i = —(3co and 
the effective index is rii = rico and the coupling constant 
K = Ki = Kco- From Eq. ^ we find that the period of 
such a grating must be comparable with the wavelength, 
Ag ^ A/2. The FWHM of a BG is usually very narrow 
compared with an LPG. Although a narrow FWHM gives 
a high resolution for a sensor, it also puts high demands 
on production tolerances and the optical spectrum ana- 
lyzer. Thus the minimal detectable shift in the resonant 
wavelength is set by the spectrum analyzer. 

Fluctuations in the effective index caused by imperfec- 
tions in the PCF structure or the grating period along 
the grating enters by Eq. ([6]) in the transmission spec- 
trum and can inflict a significant broadening of the res- 
onance peakSi. On the beneficial side one need only to 
couple light into one end of the PCF, since it is possible 
to measure the reflection spectrum of the grating. 

The coupled-mode equations for the Bragg grating may 
be solved analytically from Eq. Q yielding a transmis- 
sion coefficient 



m = 1 - 



sinh^ \k\LJ1 - 



cosh' 



{\k\L^ 



1 



52 



(11) 



At resonance the transmission is equal to T(0) = 1 — 
tanh^(|K|L). Assuming that the coupling constant is con- 
stant, K = k(Ai.), we find numerically that the ^fwhm can 
be approximated reasonably well by 



3FWHM 



4.0|k|L + 3.4 



(12) 



According to this estimate the FWHM decreases with the 
grating length until a lower bound set by the coupling 
strength itself. 

Typically \k\L > 1 and we can neglect the second term 
in Eq. and obtain (5fwhm — 4|k|. The wavelength 
FWHM is then 
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(13) 



We note that the FWHM generally increases with the 
resonance wavelength squared. 

In conclusion we have found an expression, character- 
istic for BGs, that relates Afwhm directly with the cou- 
pling constant, the group index of the core mode and 
resonant wavelength. We have made the assumption the 
coupling constant is independent of wavelength. 



4. Linear response theory for fiber sensors 

The fiber sensor is perturbed by a measurand to produce 
a change in the measured quantity, in our case the res- 
onant wavelength. For quantitative and qualitative esti- 
mation of the measurand a linear relation to the resonant 
wavelength is assumed. Since the resonant wavelength 
shifts are generally small compared to the wavelength, 
it is natural to employ perturbation theory and Taylor 
expand the resonance condition to first order around the 
resonance wavelength. 

Consider a general measurand, a. We assume that the 
shift in resonant wavelength as well as in the measurand 
are sufficiently small, and we expand the resonance con- 
dition Eq. ^ in the resonant wavelength and the mea- 
surand to first order 



A,. -I- AAr = 
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A^, we obtain an expression for the sensitivity 
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(15) 



The shift in the resonant wavelength is thus dependent 
on the change in the difference of the effective indices rel- 
ative to the difference in group indices. We have chosen 
the symbol '7' for the sensitivity to match the notation 
of Shu et al.-®. Obviously, when I7I is large the sen- 
sor will give large wavelength shifts, which will give a 
high sensitivity and accuracy of the measurand. How- 
ever, changing the fiber parameters not only changes the 
wavelength shifts but also change the full- width half max- 
imum (FWHM) of the resonance dip. Thus, even though 
I7I may be large for a certain fiber structure, the wave- 
length shifts will be difficult to detect if they are much 
smaller than the FWHM of the resonance dip. In partic- 
ular, we therefore note that the sensitivity is high close to 
group index matching, fig = 0, which by Eqs. (jlOllSp is 
also where the FWHM is wide. We introduce the quality 
factor, Q, of the grating sensor, which is the wavelength 
shift divided by the FWHM, 



Q 



1 dA, 
Afwhm da 

{ 1 L_ O(nco-nei) 
0.80 A, da ' 
^1 1 3nco 
0.64 A,|k| da ' 



LPG, 
BG. 



(16) 



We note that Q, in contrast to 7, does not increase when 
fig — > 0. Perfect group index matching fig = only 
occurs for LPGs, as we will see in the following. 

We have derived two expressions for the performance 
of a fiber grating sensor. The sensitivity, 7, is inversely 
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proportional to the group index mismatch, rig, and ex- 
presses the shift in the resonant wavelength relative to the 
resonant wavelength. The quality factor, Q, in contrast, 
expresses the wavelength shift relative to the FWHM. 
Ideally both the sensitivity, 7, and the quality factor, Q, 
should be large. 

The FWHM of LPGs is usually larger than the res- 
olution of the spectrum analyzer. The quality factor is 
thus the most important measure of the performance of 
LPGs. The FWHM of BGs is usually narrow, and com- 
parable to the resolution of the spectrum analyzer. The 
sensitivity is then the most important for BGs. 

5. Numerical simulations 

A. Effective indices 

We employ a commercial finite element code (Comsol 
Multiphysic92S) for the calculation of the effective in- 
dices. The code uses an adaptive mesh algorithm for gen- 
erating the mesh. The meshes consist of 60,000-250,000 
quadratical Lagrange elements. The core mode is simu- 
lated by considering a cross section with 12-25 rings of 
air holes: 12 rings for d/A = 0.80 — 0.95, 15 rings for 
d/A = 0.50 - 0.75, and 25 rings for d/A = 0.15 - 0.45. 
Outside the rings of air holes there is a layer of air 0.5A 
thick to reduce the effects of a finite cladding. The large 
number of rings of air holes is chosen such that effects 
of a finite cladding is small. Otherwise reflections may 
arise from both the silica-air interface and the metallic 
boundary outside the air layer. This induces an error 
on the effective index. In our simulation the number of 
rings was sufficiently large that the relative error on the 
effective indices is smaller than 10~*. Large holes give a 
strong confinement of the core mode, and obviously fewer 
rings are needed compared to small holes. 

We generally use only data where Uco — "fsm > 10^^ 
such that the error is negligible. This also makes physical 
sense, since the PCF is only weakly guiding below this 
limit. 

For the fundamental space filling mode we use the 
Bloch theorem and we need only to simulate a single unit 
cell containing only one hole. The unit cell is spanned by 
the lattice vectors seen in Fig. ([T]) of an infinite trian- 
gular cladding structure. We impose periodic boundary 
conditions on the unit cell. 

In our simulations we solve Eq. ([T]) for the propaga- 
tion constants, /?, for the modes of the fiber at a fixed 
frequency, u, which enables us to include material disper- 
sion in a consistent manner. For the refractive indices of 
the materials we use a Sellmeier expression for the silica 
and use the refractive index of water given in Refi^. All 
plots in the following are shown as function of the wave- 
length relative to the pitch, but at a fixed wavelength. 
If the material dispersion is neglected then space and 
time in Maxwell's equation may be scaled with the pitch 
of the PCF. However, in our study the material disper- 
sion cannot be neglected. We thus choose a wavelength 
A = 1.0/im. 



The grating period and therefore the beat length for 
LPGs can be determined from the effective indices us- 
ing Eq. ([6]) . The beat length varies significantly with the 
hole size and pitch, as shown in Fig. [H The resonance 
condition, Eq. ([6]), can be satisfied over a range of wave- 
lengths if OxAq = 0. For these values the FWHM of an 
LPG is very large. 

The between the curves of in Fig. |4]the gray area in- 
dicates fig < and the white area fig > 0. The boundary 
between these two areas indicates fi^ = 0. 

When the slopes of the curves in Fig. 2] are parallel 
with the normalized wavelength axis the resonance con- 
dition Eq. [6] is nearly matched for multiple normalized 
wavelengths. A/ A, we may thus expect that the wave- 
length FWHM is wide. By Eq. the FWHM is widest 
when fig = 0, indicated by the boundary of the gray area. 
It may seem surprising that the fig = does not coincide 
with zero slope of the curves. The curves in Fig. UJ how- 
ever, are evaluated for a fixed wavelength with a variable 
pitch, and not for a fixed pitch with variable wavelength. 
The difference in normalized wavelength for group index 
matching and for zero slope of the curves is caused by 
material dispersion of silica. 

For BGs fif in the resonance condition Eq. ^ only 
varies a few percent for all wavelengths, pitch values, and 
hole sizes. Thus the dispersion effects on the resonance 
condition are small. Correspondingly, it can shown that 
dispersion effects on the wavelength FWHM is also small, 
because the group index mismatch, fig, in Eq. ([8]) only 
varies slightly (< 12%) for index guiding PCFs, since Uco 
always dominates over its derivative in the group index, 

^g,co — ^co At^^TT-co- 

There is a fundamental difference between PCFs and 
standard optical fibers is that a PCF consists of two, of- 
ten very different, materials. Standard optical fibers also 
consist of two materials: pure silica and a core of doped 
silica. In practice, however, the material parameters for 
the doped silica may often, to a good approximation, be 
taken to be equal to those of pure silica. For a PCF, 
e.g. the strain optic response of silica is completely dif- 
ferent from that of air, which is obviously zero for our 
applications. 

To correctly handle a two material fiber one must de- 
termine the field energy intensity fraction, /„, inside the 
holes of the PCF. The field energy fraction for the base 
material is then, 1 — /«. It is not possible to simply add a 
contribution, e.g. stress optic, directly to the resonance 
condition, Eq. ([6]). Instead every contribution must be 
weighted with the fraction of field intensity in the corre- 
sponding material, as shown in the appendix. 

It is possible to make some generalizations about field 
energy fraction in the holes, /„, for index guiding fibers. 
In the short wavelength limit where the wavelength is 
much smaller than the pitch of the structure, A <C A, the 
field is confined to the high index base material. In the 
long wavelength limit, A ^ A, the changes in the field 
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are small over a length of the pitch. This imply 



6. Refractive index sensing and biosensing 



/„^0 for A/A^O, 
fu^ f for A/A -> oo, 



(17) 
(18) 



where / is the air filling fraction, i.e. the fraction of area 
covered by the holes in the cross section of the fiber. The 
approximation /« = is natural as a first approximation 
since the field intensity in the air holes is generally 
much smaller than the field intensity in the base material. 
But PCFs can possibly have air filling fractions of almost 
unity, implying that also /„ can be almost unity, and 
therefore the effect can not be neglected. 

We have shown /„ for the FSM in Fig. [S] for air filled 
holes {rih = 1.0). /„ for the core mode of a PCF with 
water filled holes {uh — 1.33) is shown in Fig. ??) and 
the corresponding FSM is shown in Fig. [51 Clearly /„ 
is larger for the FSM than for the core mode, and it 
increases with increasing wavelength and increasing hole 
diameters d/A. The /„ also increases when the index 
contrast between the hole and base material is decreased. 
Hence is larger for water filled PCFs than for air filled. 

B. Small and large core PCF 

In the following we refer to two PCFs to as the large core 
and the small core PCF. The structure parameters for the 
small core PCF is (d/A, A) = (0.45, 7^m) 2) and for the 
large core PCF (d/A, A) = {0.7b, 2fim). We choose the 
length of the gratings to be 30 mm for LPGs with addi- 
tionally K — 4/L for BGs. The sensitivity, 7, is indepen- 
dent of these parameters. The quality factor, Q, depends 
inversely proportional on the length of the grating, when 
assuming constant kL. For BGs the quality factor also 
depends inversely proportional to the coupling constant, 



A. Refractive index sensing 

Refractive index sensors can be used in the application 
of continuous monitoring of gasses in the oil industry. 
The gas flows through the PCF is detected by measuring 
the overall refractive index of the hole contents. Opti- 
cal fibers have the advantage that there are no electric 
currents thus avoiding completely the possibility of gen- 
erating electric sparks. 

The sensitivity for refractive index sensing is also 
used for benchmarking the sensitivity of evanescent wave 
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Fig. 7. Field energy intensity fraction in the water filled 
holes for the fundamental space filling mode. The lines 
indicate different values of the hole diameter relative to 
the pitch: 0.15 to 0.95 in steps of 0.10 increasing from 
below 
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Fig. 6. Field energy intensity fraction in the air filled 
holes for the fundamental space filling mode. The lines 
indicate different values of the hole diameter relative to 
the pitch: 0.15 to 0.95 in steps of 0.10 



Fig. 8. Field energy intensity fraction in the water filled 
holes for the core mode. The lines indicate different val- 
ues of the hole diameter relative to the pitch: 0.15 to 0.95 
in steps of 0.10 increasing from below 
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biosensors, which will be studied in the next section. 

Alternatively, a refractive index sensor may be used as 
tunable grating if the holes are filled with a substance 
whose refractive index can be controUed^'''^^'^''. Tunable 
gratings may have applications within modulation of op- 
tical signals. 

Infiltrating substances into PCF-LPGs that are ab- 
sorbing at the resonant wavelength requires special con- 
siderations, as discussed by Daxhelet et al.— . In brief the 
attenuation induced by the absorption of the substance 
becomes significant if the modes are significantly atten- 
uated over the length of a period of the grating. The 
interference of modes in the grating becomes weakened 
and the resonance dip is broadened. Since the grating 
period of an LPG is much longer than that of a BG, the 
LPG can be considerably affected by absorbtion, while 
the BG is largely unaffected even by a high absorbtion. 

The change in the effective index of a mode induced by 
a change of refractive index of the hole contents is pro- 
portional to the fraction of the time averaged field energy 
intensity {u{r±) = ^e{r ±)\Ei{r ±)\'^) inside the holes. For 
BGs the change in flf = 2nco is (see appendix for deriva- 
tion) 

|^ = 2!^/_, (19) 

where An/j is the change in the refractive index, n^, of 
the content in the hole. fu,co is the fraction of the field 
energy intensity of the core mode inside the holes. The 
sensitivity is also proportional to the group index ng,co- 
For LPGs the change in hf = rtco — 't-cI is 

o — Ju, CO Ju,c\- y/'^ } 

OUh Ufi nil 

It is clear from Eqs. (|15ll6p that the sensitivity and the 
detectability are both largely determined by /„. For BGs 
fu should be as large as possible. 

For BGs the group index is close to the effective in- 
dex, rig, CO — — Xdxrico ~ rico since the derivative is 
negligible as seen in Fig. O Thus for BGs the sensitivity 
is determined by /„,co- In Fig- [5] we realize that a large 
/u.co implies a small pitch and a large hole diameter. For 
BGs the size of Q is largely determined by the pitch. Q 
also increases with decreasing hole size, but to a lesser 
extent. 

In order to have a high Q for LPGs the difference 
%,co/«,co - '^g,ci/u,ci should be large. Comparing Figs. [8] 
and [7] we realize that this also implies a small pitch and 
a large hole diameter. But in this case the difference falls 
off beyond a certain point in the wavelength pitch ratio. 
This is because /u,co and /,j,ci approaches the same limit 
for large wavelengths, and thus the Q decreases, as seen 
in Fig. M 

In Fig. [9] there is also indicated the group index 
matched line for LPG where nf = Hg,co — JT-gjCi = 0. Along 
this line the peak of the LPG is very wide by Eq. ([5]). On 
the left side of this line the shift in wavelength are positive 



for a positive Ann- On the right hand side the shift in 
wavelength is negative for positive Auh- The reason for 
this, is that fig enters in Eq. (|15p in the denominator, and 
thus its sign determines the sign of the wavelength shift. 
For BGs group index matching nf = 2nco — 2A9a?t.co = 
does not occur since for the considered index guiding 
PCFs the effective index is always larger than its deriva- 
tive. 

For the sensitivity of BGs Eq. (|19p is inserted in 
Eq. P3|) . and plotted for the wavelength 1.0 iim, in 
Fig. [TUl For LPGs Eq. (PO)) is inserted in the expres- 
sion for the Q, Eq. (fT6|) and plotted for the wavelength 
1.0 Aim, Fig.d 

To give some specific details we choose a typical large 
mode area, endlessly single mode PCF with structure 
values: A ~ 7^m and d/A ~ 0.45. This gives a core 
diameter of 10 /im. We choose the length of the grating 
to be 30 mm, and the coupling strength is k = tt/{2L). 
We choose different wavelengths for both LPGs and BGs: 
600, 900, and 1550 nm. The data is presented in Table 
[T] The sensitivity increases with increasing wavelength 
since the evanescent wave is increasing inside the holes, 
corresponding to an increasing /„. The sensitivity, 7, for 
LPGs is four orders of magnitude larger than the sensi- 
tivity for BGs. But taking the FWHM of the resonance 
peak into account, the Q is only slightly larger for LPGs 
than for BGs. 

The data for a small core PCF are presented in Table 
[T] The small core PCF has significantly higher sensitivity 
than the large core PCF. Notice that the sensitivity for 
the small core has opposite sign to the large core PCF. 
Also Q decreases with increasing wavelength. The shift is 
negative because the PCF's structure parameters places 
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Fig. 9. Quality factor Qri for LPG refractive index sens- 
ing. The lines indicate different values of the hole diam- 
eter relative to the pitch: 0.25 to 0.95 in steps of 0.10. 
The gray area indicates negative group index mismatch 
ng < 0. The wavelength is 1 ^m, the length of the grat- 
ings is L = 30 mm, and Ik^I = 7r/(2L) 
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large-core PCF 


Refractive index sensing 


Biosensing 


-^FWRM 


7RI 


Qri 


AA/Arih [nm] 


7Bio [1/nm] 


Qbio [1/nm] 


AA/fBio 


LPG, 600 nm 


1.6 


43 


980 nm 


2.4x10^ 


63x10"^ 


1.4 nm/nm 


23 nm 


LPG, 900 nm 


3.2 


88 


2700 nm 


3.0 xlO^ 


85x10"^ 


2.7 nm/nm 


32 nm 


LPG, 1550 nm 


14 


210 


22x10^ nm 


8.1 xlO"^ 


120x10"=* 


12 nm/nm 


109 nm 


BG, 600 nm 


0.10x10"^ 


5.8 


0.060 nm 


3.0x10"'^ 


17x10"^ 


0.18 pm/nm 


10 pm 


BG, 900 nm 


0.34x10"^ 


13 


0.30 nm 


6.5x10"'^ 


25x10"^ 


0.59 pm/nm 


24 pm 


BG, 1550 nm 


1.6x10"^ 


35 


2.5 nm 


19x10"^ 


41x10"^ 


2.8 pm/nm 


70 pm 


small-core PCF 


Refractive index sensing 


Biosensing 


^FWHM 


7ri 


Qri 




7Bio [1/nm] 


Qeio [1/nm] 


AA/fBio 


LPG, 600 nm 


73 


6500 


44x10^ nm 


89x10"^ 


7.9 


-54 nm/nm 


6.8 nm 


LPG, 900 nm 


11 


8600 


9.8x10^ nm 


7.7 xlO"^ 


6.2 


-7.0 nm/nm 


1.1 nm 


LPG, 1550 nm 


6.4 


7400 


9.9x10^ nm 


2.7 xlO~^ 


3.0 


-4.1 nm/nm 


1.3 nm 


BG, 600 nm 


6.0x10^^ 


350 


3.6 nm 


17x10"^ 


1.0 


10 pm/nm 


10 pm 


BG, 900 nm 


18x10"^ 


690 


16 nm 


35x10-*^ 


1.33 


31 pm/nm 


23 pm 


BG, 1550 nm 


74x10"^ 


1700 


115 nm 


83x10"^ 


1.9 


130 pm/nm 


70 pm 



Table 1. Comparison of sensitivity and the quality factor for PCF biosensors and refractive index sensors for both 
LPGs and BGs at different wavelengths for two different PCFs: a large core of 10 /im and a small core PCF of 1.5 
/xm. 



the PCF on the left of the group index matching line 
ng = indicated in Fig. [9l 

The minimum detectable change in refractive index 
for LPGs can be calculated. For L = 30 mm, we get 
a maximal Q ^ 10^. Taking that the minimal detectable 
wavelength is ApwHivi/lOO we get a maximal sensitivity 
of Anhanin = j^q ^ 10"^. The same calculation for a 
BG with Q ^ 10^ gives Anh,min ~ lO""*. 

Both LPGs'^" and BGs^^ have been used for measuring 
the refractive index of a liquid analyte inside the holes 
of the PCF. Comparing with experiment for LPGs2£ the 
PCF has parameters d/A = 0.47, A = 7.2/ini, and the 
LPG has Ar — 845nm and L = 18.2 mm. The minimal 
detectable wavelength shift is AAmin ~ 1 nm. Theo- 
retically the parameters give Qm — 52 (no units) and 
Afwhm — 45 nm which agrees with the experimental 
value ApwHM ~ 45 nm. The minimal sensitivity given 
the mentioned minimal detectable wavelength is theoret- 
ically Anh,min ~ 4.3 x 10"'^ which is in agreement with 
the experimental value Anh^min ~ 10"^. The shift m 
resonances wavelength as function of sensitivity is exper- 
imentally AAr/Arih — 7500 nm but gives theoretically a 
lower value AAr/Arth — 2300 nm. 
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Fig. 10. 7RI for BG refractive index sensing. The lines 
indicate different values of the hole diameter relative to 
the pitch: 0.15 to 0.95 in steps of 0.10. The wavelength 
is 1 /im, the length of the gratings is L = 30 mm, and 
\kL\ = 4 



For BGs we compare with Huy et ali^. They have a 
6 hole fiber with c? « 15 /im and pitch « 15.8 /im. This 
gives theoretically Afwhm — 117 pm in agreement with 
the experimental Afwhm ~ 150 pm. Theoretically we 
have AAr/Anh ~ 0.59 nm and taking their a minimal 



detectable wavelength shift of AAmin = 1-0 pm, we find 
Anh,min ~ 1-7 x 10""^ which is in agreement with the 
experimental value Anh,min ~ 4 x 10~^. Their numerical 
simulations gives a minimal detectable shift which is a 
magnitude off the experimental value. 
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B. Label-free biosensor 

Biosensors in PCFs have seen some interesting devel- 
opments in recent times. The PCF biosensors can be 
made to detect both DNA, antibodies or antigens, and 
the PCF can be made of both siUca and polymer materi- 
als. The experiments have been conducted using labeled 
molecules, i.e. target biomolecules have a fluorescent tag, 
that can be detected using fluorescence or absorbance mi- 
croscopy. To be really interesting the measurements must 
be label-free. This can be done with PCF-LPG using the 
evanescent-wave sensing principle. This principle is used 
in the widespread surface-plasmon resonance biomolecule 
interaction analysis instruments. 

The evanescent-wave principle measures the increase 
of refractive index at a surface due to the presence of 
biomolecules. The experiments are typically carried out 
in an aqueous environment. Water has refractive in- 
dex ~ 1.33 and biomolecules have ^ 1.45 — 1.48. A 
given biomolecule may be tested for its affinity to an- 
other biomolecule, by immobilizing the first biomolecule 
onto the surface (functionalization), and then infiltrat- 
ing the second molecule into the PCF. By a wash any 
biomolecules that have not been immobilized may be 
flushed out. The immobilized biomolecules will consti- 
tute an increase in the layer of molecules on the surface. 
The layer is 1-10 nm thick depending on the application. 

The film of molecules on the surface may not consti- 
tute a homogeneous layer. But since the wavelength is 
much longer than the thickness of the layer (1 micron ^ 
10 nm), the sensor measures the thickness of the biofilm 
on average. The sensor is thus largely insensitive to in- 
homogeneous distributions of molecules on the surface. 
For benchmarking we assume that the layer of molecules 
has a homogeneous thickness ^Bio and has a refractive 
index of 1.45. This refractive index is to a good ap- 
proximation identical to silica: rir^Bio — 'T-Si02- Thus a 
layer of biomolecues ^Bio thick corresponds to decrease 
of Ad = — tBio/2 in the hole diameter. To calculate the 
change in effective index due to a layer of molecules we 
thus find 

dnes _ dnes , . 

a^Bio ~ dd' ^''> 

where d is the hole diameter. The change in fi{ for LPGs 
and BGs is 

dtBio \ -4^, BG. ^ ' 

In Figs. [Tl] and [12] it is seen that Qbio and 7Bio are 
largest for small pitches and large wavelengths. For BGs 
the 7Bio is largely governed by the pitch. For LPGs a 
large benefit comes from using large holes instead of small 
holes, whereas the pitch is less important. 

For BGs Eq. (|22p is inserted in the expression for the 
sensitivity, Eq. p3|) and plotted for the wavelength 1.0 
fim, Fig. [T^l For LPGs Eq. (22) is inserted in the expres- 
sion for the Q, Eq. (fTB|) and plotted for the wavelength 
1.0 fim, Fig.[nj 



To compare with experiments we take the PCF struc- 
tural values from^: d/A = 0.47, A = 7.2^m. The 
wavelength was « 845 nm. This gives AA — 2.3 x tBio, 
Q = 53 X 10-3 (no units), and 7 = 2.8 x lO^^/nm. Thus 
a 1 nm thick layer of biomolecules gives a 2.3 nm shift in 
resonant wavelength. It should be noted that the paper- 
gives a value (1.4 nm/nm) that is ~ 50% more. This is 
because the effective indices were calculated using a dif- 
ferent numerical method, and material dispersion was not 
included in the simulations. 
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Fig. 11. Quality factor, Qbio, for LPG biosensing. The 
lines indicate different values of the hole diameter relative 
to the pitch: 0.25 to 0.95 in steps of 0.10. The gray area 
indicates negative group index mismatch fig < 0. The 
wavelength is 1 fim, the length of the gratings is L 30 
mm, and \kL\ = n/(2L) 
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Fig. 12. Sensitivity, 7Bio, for BG biosensing. The lines 
indicate different values of the hole diameter relative to 
the pitch: 0.15 to 0.95 in steps of 0.10. The wavelength 
is 1 /im, the length of the gratings is L = 30 mm, and 
\kL\ = 4 
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In comparison the 6 holes PCF of Huy et ali^i (d « 
15 ^m, pitch « 15.8^m) gives AA — 0.73 x 10~^ x ieio, 
Q = 6.2 X 10"3/nm, and 7 = 0.47 x 10"*^/nm at A,. ~ 
1555nni. Thus a 1 nm thick layer of biomolecules gives a 
0.73 pm shift in resonant wavelength. 

In general, both the sensitivity and the quality fac- 
tor for BGs is monotonically increasing with increasing 
wavelength. They also increase with increasing hole size 
as seen in Table [TJ The general behavior for the both 
the sensitivity and the quality factor for LPGs is more 
complicated. The sensitivity is strongly influenced by 
the group index matching, and the resonant wavelength 
shifts can change sign depending on the PCF structure 
parameters. The quality factor, Q, is high for large air 
holes, and big normalized wavelength A/ A. 

In conclusion, the sensitivity of label-free PCF grating 
biosensors follows roughly the tendencies of the refractive 
index sensors. The sensitivity varies several orders of 
magnitude in sensitivity and quality factor with the PCF 
structure parameter. BGs have a high sensitivity when 
the pitch is small compared with the wavelength. Thus 
the optimal BG biosensor sensor is a small core fiber. 
For LPGs the quality factor is high when the air holes 
are large, and are less influenced by wavelength to pitch 
ratio, except in the short limit where the quality goes 
to zero. The optimal LPG biosensor or refractive index 
sensor thus has very large holes and a pitch about three 
times or less the wavelength. 



7. Temperature and strain sensing 

A. Temperature 

Sensitivity to temperature is often an unwanted source 
of error in fiber gratings. However, fiber gratings may 
also be employed as high temperature sensors since silica 
is stable to high temperatures. If the susceptibility to 
temperature is an unwanted effect it may also be reme- 
died by using two multiplexed fiber gratings each with 
a different linear temperature response. In the following 
we will study the temperature response of air filled sil- 
ica PCFs although the theory is also applicable to other 
materials. 

Temperature changes both the material refractive in- 
dices of the PCF through the thermo-optic effect and 
also by expansion of the PCF structure. The thermo- 
optic coefficient for the refractive indexof silica glass is 
~ 9 X 10^®/K and the thermal expansion coefficient 
is a ~ 0.55 X 10^^/K. Since the perturbation of the PCF 
is sufficiently small we may expand the effective index 
to first order in the temperature. Since both the pitch, 
length and the refractive index of the silica are dependent 
on silica, we use the partial derivatives of the effective re- 



fractive index (the chain rule) 



dn^s 



dT 



dneff(i(T),A(T),n,(r)) 



drics dL 



dT 

ducff OA 



dL dT dA dT dnu dT 



(23) 



We readily identify the thermal expansion^ ~ aA and 
^ = aL. The partial derivative of the effective index 
with respect to the length can be evaluated according to 
following argumentation. An elongated PCF at a fixed 
effective index is equivalent to a PCF at fixed length with 
an increased effective index. Formulated in mathematical 
terms this yields 



dL 



(24) 



In Eq. (|23p the partial derivative of the effective index 
with respect to the pitch (2nd term rhs) must be eval- 
uated numerically. The contribution from the thermo- 
optic effect can be evaluated using Eq. Eq. (1^ 
can then be written for a BG as 



dn[ 
'df 



2a(nr 



A 



X, 



~dA 



For an LPG we get 



dui 
'df 



— Q!(nf-f A 



driis 
'dA' 



+- 



b,T 



rib 



[(1 - fu,co)ng 



/u,co)^g,cc(-25) 

(26) 

(1 - /«,cl)ng,cl] ■ 



For generality we have kept a general refractive index of 
the base material. 

For BGs Eq. ((25|) is inserted in the expression for the 
sensitivity, Eq. (fT5|l . and plotted for the wavelength 1.0 
/im. Fig. [Ml For LPGs Eq. ([26]) is inserted in the expres- 
sion for the Q, Eq. (fTB)) . and plotted for the wavelength 
1.0 ^lm in Fig.[ll 

The sensitivity of BG PCF temperature sensors is al- 
most constant for the structure parameters, A and d/A, 
Fig. HH The Q value for PCF-LPG temperature sensors 
in contrast varies several orders of magnitude with the 
structure parameters, A and d/A, Fig. [T31 The dips 
seen on the left in the figures is a change of sign of the 
wavelength shift. 

In table [5] we see that the resonant wavelength shifts 
for both BGs and LPGs are of the order of 5pm/K, re- 
gardless of PCF parameters. The quality factor is sev- 
eral orders of magnitude smaller for LPGs compared with 
BGs. This can be attributed to the a partial cancelation 
of the thermo-optic terms of the core and cladding mode 
in Eq. 

We note that the sensitivities for the BGs are surpris- 
ingly constant in Table d In the rhs of Eqs. (HI [HD 

the terms ricff + A^ and ng(l - /„)/nf, are roughly 
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of the same order of magnitude. The prefactor for the 
thermo-optic term is larger than the prefactoc of ther- 
mal expansion term, Cj/a w 20. The temperature re- 
sponse of the grating is largely determined by the thermo- 
optic. For BGs we may thus approximate 7^ = l/X^f- — 
Clln^,I\T. For silica PCFs 7t ^ Cljn^, ~ 0.62x 10.6/K 
in good agreement with the values given in tables [2] and [21 
We conclude that sensitivity of BG temperature sensors 
is constant, within a few % margin, for all PCF struc- 
tures and is dominated by the thermo-optic cocfhcicnt 
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Fig. 13. Quality factor, Qt, for LPG temperature sens- 
ing. The lines indicate different values of the hole diam- 
eter relative to the pitch: 0.25 to 0.95 in steps of 0.10. 
The gray area indicates negative group index mismatch 
fig < 0. The wavelength is 1 /im, the length of the grat- 
ings is L = 30 mm, and \kL\ = ■k/(2L) 
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relative to the refractive index of the base material. 



To validate our theory we compare with experiments. 
Dobb et al.'^^ characterized LPGs in PCFs inscribed with 
an electric arc. For a PCF with structure parameters 
A ~ 7.1 /xm, d ~ 4.0 ± 0.4 ^m they found a shift of 
resonance wavelength with temperature as 3.4 pm/K at 
the wavelength 1511 nm. Using the structure parameters 
in the theory yields 4.3 pm/K. For a PCF with structure 
parameters A ~ 8 /im, d ~ 3.7/im, they found 2.2 pm/K 
at 1403 nm wavelength and the theory gives 6.1 pm/K. 
Thus we have an agreement of the right magnitude and 
right sign of the wavelength shift. It is seen in Fig. [T3] 
that for their structure values, d/A ~ 0.46, A/A ~ 0.18 
the quality factor, Q, experiences a series of dips. The 
dips correspond to the terms in Eq. (pS)) cancel in the net 
contribution. Both numerically, and experimentally, an 
estimation of the sensitivity will be prone to small errors 
on contribution from each term. 



In an earlier experiment Humbert et dl A used the elec- 
tric arc method to inscribe LPGs in a PCF with structure 
parameters A ~ 4.0/im, d ~ 2.0/im. They found a value 
of 9 pm/K at 1590 nm, roughly half the theoretical value 
-23 pm/K, but with an opposite sign. To inscribe BGs 
Martelli et al. used a PCF with a erbium-doped UV 
sensitive core^. Their structure parameters are roughly 
A ~ 10/im, d ~ 2/xm. At a wavelength of 1535 nm the 
wavelength shift is 20 pm/K. The theory gives 10 pm/K 
for a pure silica PCF. An erbium doped core is known to 
decrease the temperature sensitivity in standard optical 
fibers. 



BGs have also been inscribed in a small core fiber with 
the structure parameters d/K = 0.5 and A = 1.6/im^^. 
The core was doped with germanium to make the core 
UV-sensitive. They found an experimental sensitivity of 
4.11 pm/K at a wavelength of 1508 nm. Our theory 
gives 9.9 pm/K. A Germanium doped core is known to 
decrease the temperature sensitivity in standard optical 
fibers. 



Fig. 14. Sensitivity, 7t, for BG temperature sensing. 
The lines indicate different values of the hole diameter 
relative to the pitch: 0.15 to 0.95 in steps of 0.10. The 
wavelength is 1 /im, the length of the gratings is L = 30 
mm, and \k.L\ ~ 4 



To enhance the temperature sensitivity in PCF the 
holes can be filled with substances with a high thermo- 
optic coefficient. This has been shown to increase tem- 
perature sensitivity^^. Eq. ((26)) can in this case be modi- 
fied to be valid by adding an additional temperature optic 
term for the holes. 
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large-core PCF 


Temperature 


Strain 


-^FWHM 


7t 


Qt [1/K] 


A A/ AT [1/K] 






AX/es 


LPG, 600 nm 


3.9x10"® 


UxlO"'"' 


2.3 pm 


-1.2x10"" 


4.5 xlO"-"^ 


-0.74 pm 


17 nm 


LPG, 900 nm 


6.2x10"" 


30x10"^ 


5.6 pm 


-2.3x10"" 


11x10"^ 


-2.1 pm 


19 nm 


LPG, 1550 nm 


1.7x10"" 


10x10"^ 


27 pm 


-2.1x10"" 


13x10"^ 


-3.3 pm 


26 nm 


BG, 600 nm 


6.7x10"" 


0.39 


4.0 pm 


0.82x10"" 


0.047 


0.49 pm 


10 pm 


BG, 900 nm 


6.8x10"" 


0.26 


6.1 pm 


0.82x10"" 


0.031 


0.73 pm 


24 pm 


BG, 1550 nm 


6.8x10"" 


0.15 


11 pm 


0.81x10"" 


0.018 


1.3 pm 


70 pm 


small-core PCF 


Temperature 


Strain 


-^FWHM 


7t 


Qt [1/K] 


AX/ AT [1/K] 






AA/£, 


LPG, 600 nm 


5.0x10"^^ 


5.8x10"^ 


-3.0 pm 


-2.4x10"" 


2.8x10"^ 


-1.4 pm 


0.52 nm 


LPG, 900 nm 


4.0x10"^ 


23x10"^ 


-36 pm 


-4.2x10"" 


2.4x10"^ 


-2.0 pm 


1.6 nm 


LPG, 1550 nm 


3.7x10"^ 


55x10"^ 


57 pm 


0.7x10"" 


1.0x10"^ 


1.1 pm 


1.0 nm 


BG, 600 nm 


6.7x10"" 


0.39 


4.0 pm 


0.81x10"" 


0.047 


0.48 pm 


10 pm 


BG, 900 nm 


6.7x10"" 


0.26 


6.1 pm 


0.79x10"" 


0.030 


0.71 pm 


23 pm 


BG, 1550 nm 


6.6x10"" 


0.15 


10 pm 


0.75x10"" 


0.017 


1.2 pm 


68 pm 



Table 2. Comparison of sensitivity and the quality factor for PCF temperature and strain sensors for both LPGs and 
BGs at different wavelengths for two different PCFs: a large core of 10 /im and a small core PCF of 1.5 /im 



B. Strain 

Fiber optic strain sensors have found application within 
civil engineering, where they provide a stable and sensi- 
tive measurement of the strain in constructions. Strain 
is defined as the relative elongation: Eg = AL/L of a 
piece PCF of length L. Often millistrain, fiSg = lO^e, is 
used instead of strain. When the PCF is elongated, the 
PCF contracts transversely to the applied strain. This is 
expressed by Poisson's ratio which formally is the ratio 
between lateral and tensile strain with a negative sign. 
For silica glass Poisson's ratio is given by i/ = 0.17. Strain 
also changes the refractive index of materials. The strain- 
optic coefficient for silica is = rj^ — —0.26^. Since 
strain changes the length, pitch, and the material refrac- 
tive indices we apply partial derivatives (the chain rule) 
to the linearization of the effective index 



dn^a dL 
dL dss 



de^ 

ducff dA 



dA dss dur dss 



(27) 



Since L{es) 



(1 -I- es)L we have ^ 



L and A{£s) 



(1 + £s)A gives ^ — —vA. An increase in length of 
the PCF with a fixed effective index is mathematically 
equivalent to a fixed length with an increase in effective 
index. 



dL 



"cff 

L ' 



(28) 



The derivative of the effective index with respect to the 
refractive index of the base material may be evaluated by 
inserting the strain-optic coefficient in the perturbative 
formula, Eq. (|36p. given in the appendix. The formula 
for BGs is 



2nr 



vA 



d2nc( 
dA 



2^g,co(l /tt,co), 



while for LPGs we obtain 



{ric 



f^cl) 



vA 



d{nc 



nci) 



dA 



(29) 



(30) 



+ — ("-g,co(l - /u,co) - "-g,cl(l - /u,cl))- 
Tib 



Eqs. (|29l30p are inserted in the expression for the sensi- 
tivity, Eq. (fT5|) and calculated for the wavelength 1.0 /im. 
The results are shown in Figs. [16] and ??, respectively. 
The sensitivity for BGs is very close to being constant re- 
gardless of pitch, hole diameters, and wavelength, as seen 
in Fig. [TBI and Table [21 This can understood by that fact 
that Hco and ?T.g,co are of similar magnitudes, such that 
TT-co ^ A^^. We can neglect the second term in Eq.[^ 
We then put Anf/es = (1 -I- %/nh)ng^co — 0.8 x Wg^co, 
where rtf, ~ 1.45 is the refractive index of silica. Inserting 
this in the expression for the sensitivity, Eq. ([T5|), we real- 
ize that 7e,//ies « 0.8 X 10"^ and AA = 0.8 x 10"^ A/xe^. 
This result is in very good agreement with the values for 
BGs displayed in table [2] and in Fig. [TBI 
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The quality factor, Q^^, for LPGs with A/A > 0.2 
tends to increase for increasing hole size, while being al- 
most independent of wavelength, as seen in Fig. [151 For 
these wavelength the first term in Eq. the elon- 

gation of the PCF, is dominating over the strain-optic 
effect and the lateral strain. For A/A < 0.2 the sensi- 
tivity is generally lower since fif decrease for decreasing 
wavelengths, but also because the strain-optic effect and 
the lateral strain become significant. The dips indicate a 
net cancelation of the terms in Eq. (|30p . 



10" 




10 I • • • • 1 

0.2 0.4 0.6 0.8 1 

X/A 



Fig. 15. Quality factor, , for LPG strain sensing. The 
lines indicate different values of the hole diameter relative 
to the pitch: 0.25 to 0.95 in steps of 0.10. The gray area 
indicates negative group index mismatch fig < 0. The 
wavelength is 1 fim, the length of the gratings is L — 30 
mm, and [k^I = tt/{2L) 







8.8 




8.6 


d/A = 0.25 


8.4 




8.2 




8 




7.8 




7.6 




7.4 


d/A = 0.95 


7.2 




■7 





0.2 0.4 0.6 0.8 1 

X/A 



Fig. 16. Sensitivity ,7e^, for BG strain sensing. The lines 
indicate different values of the hole diameter relative to 
the pitch: 0.15 to 0.95 in steps of 0.10. The wavelength 
is 1 jim, the length of the gratings is L = 30 mm, and 
|kL| =4 



To validate our theory we compare with experiments. 
Dobb et al,''^ characterized LPGs in PCFs inscribed with 
an electric arc. For a PCF with structure parameters 
A ~ 7.1/im, d ~ 4.0/im they found a shift of reso- 
nance wavelength with -2.5 pm//i£s at the wavelength 
1688 nm. Using the structure parameters in the theory 
yields -3.7 pm/fiSs- For a PCF with structure param- 
eters A ~ 8fim, d ~ 3.7/im, they found -2.08/ fiEs at 
1403 nm wavelength and the theory gives -3.2/iies- Thus 
we have an agreement of the right magnitude and right 
sign of the wavelength shift. Their PCF structure values, 
d/A ~ 0.46, A/A ~ 0.18, gives a Q that is close to the dips 
in Fig. [151 The dips corresponds to a cancelation of the 
terms in Eq. ([50)) and close to the dips the quality factor 
is more susceptible to errors between experimental and 
numerical parameters. Martelli et al. used a PCF with a 
erbium-doped UV sensitive core to inscribe BGs^. Their 
structure parameters are roughly A ~ 10/im, d ~ 2/im. 
At a wavelength of 1535 nm the wavelength shift is 1.2 
pm/fiSs- The theory gives 1.25 pm/^e^ for a pure sil- 
ica PCF. A doped core does not change the strain-optic 
coefficient significantly due to a doped core^. 

For BGs in a small core fiber with the structure pa- 
rameters d/A ~ 0.5 and A ~ 1.6/imi^ it is found exper- 
imentally a sensitivity of 1.25 pm//Lt£ at the wavelength 
1508 nm. The core was doped with germanium to make 
the core UV-sensitive. Our theory gives 1.16 pm//ie in 
good agreement with the experimental value. 

In conclusion we have found that the sensitivity for 
BGs is approximately constant both for temperature and 
strain sensing, with the coefficients jt — 0.62 x 10.6/K 
and 7^e3 — 0.8 x 10~^ for pure silica PCFs. The wave- 
length shift both for BGs and LPGs is generally of the 
order of 5 pm/K, and the sensitivity and quality fac- 
tor are an order of magnitude large (in K) compared 
with strain (in microstrain). The quality factor for LPGs 
for temperature and strain sensing has roughly the same 
characteristics. 

8. Conclusion 

In conclusion we have shown that the performance of 
LPG sensors are best characterized by the detectabil- 
ity, Q, rather than the sensitivity, 7. We have calcu- 
lated expressions for the sensitivity, that incorporates the 
material-optic responses to the external perturbation, as 
well as treated effects of having air holes in the struc- 
ture. Using the results of a fuUy-vectorial finite element 
calculations, incorporating material dispersion, the sensi- 
tivity and detectability have been calculated for a range 
of PCF structure parameters, i.e. hole diameters and 
pitch lengths. 

We have proposed expressions for the sensitivity for 
internal refractive index sensing, label-free biosensing, 
strain, and temperature sensing. For strain and tem- 
perature we have included both thermo-optic and strain- 
optic effects, as wells as the effects of both longitudinal 
and lateral elongation. The formulas are of general form 
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and apply to all PCFs as well as standard optical optical 
fibers. 

For PCFs infiltrated with samples, the Q for refrac- 
tive index sensing and biosensing may be improved by 
several orders of magnitude by appropriate choice of the 
hole size and pitch. For temperature and strain index 
sensing with BG sensors we have shown that the sen- 
sitivity is almost constant for all PCFs. For LPGs an 
accurate theoretical prediction of the sensitivity for tem- 
perature sensing or strain sensing can be difficult since 
the sensitivity strongly affected by the waveguide disper- 
sion. However, our analysis shows that the low tempera- 
ture sensitivity of PCF-gratings, ~5 pm/K, is generally 
valid for all PCFs. 

L. Rindorf 's e-mail address is lhr@com.dtu.dk. 

9. Appendix: First order perturbation theory 

The linear change in propagation constant by a small 
perturbation in the dielectric function, Ae, is^ 



A/3„ 



ojd(3 /^rfr^A£(r^)|E„,(r^)|^ 
2duj ^drie(ri)|E„,(ri)|2 



(31) 



Q. denotes the fiber cross section. By setting the norm of 
the orthogonal solutions to the Helmholtz equation ([1]) 
to unity. 



/ dr^Ej„(r_L)e(r_L)E„(rj^) = (5^„, 



(32) 



where f denotes Hermitian conjugate (the complex con- 
jugate transposed), we get 



Ant 



extrmef f ,m 



^^dr^Ae(r^)|E™(r^)|2 (33) 



with the group index, n 



We define the field 



energy intensity in the air holes of the mode m by 



fu 



/^driDL(ri) -E 



(34) 



is the part of the fiber cross section covering holes, 
and D — eE. The linear change in effective index of 
mode m caused by a change of the refractive index of the 
holes is then 



9nh 



fu,i 



(35) 



where nh is the refractive index of the holes. The change 
is proportional to fu,m The corresponding linear change 
in effective index caused by a change of the refractive 
index of the base material is 



due 



duh 



"(1 /u,m) 



(36) 



where is the refractive index of the base material. It is 
seen that the change is proportional to the field intensity 
in the base material, 1 — fu- 
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